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ABSTRACT. In this paper we are concerned with finding the vertices of the 
Voronoi cell of a Euclidean lattice. Given a basis of a lattice, we prove that 
computing the number of vertices is a # P-hard problem. On the other hand we 
describe an algorithm for this problem which is especially suited for low dimen- 
sional (say dimensions at most 12) and for highly-symmetric lattices. We use 
our implementation, which drastically outperforms those of current computer al- 
gebra systems, to find the vertices of Voronoi cells and quantizer constants of 
some prominent lattices. 



Let L = BJ/ 11 C W 1 be a lattice of rank m in Euclidean space given by a 
matrix B G R nxm of rank m. By linL we denote the linear subspace spanned by 
the elements of L. The Voronoi cell of L is 



The Voronoi cell of a lattice is a centrally symmetric, convex polytope. The poly- 
topes V(L) + v for v G L tile lin L. The study of Voronoi cells is motivated by the 
fact that most important geometric lattice parameters have a direct interpretation in 
terms of the Voronoi cell: The determinant det L equals the volume of V(L), the 
packing radius \{L) equals the inradius of V(L), the covering radius fi(L) equals 
the circumradius of V(L), and the quantizer constant G(L) is 



In this paper we consider theoretical and practical aspects of the computation 
of the covering radius as well as the quantizer constant of a lattice. These two pa- 
rameters have many applications, we just name a few: By computing the covering 
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1. INTRODUCTION 



V(-L) = {x G lin L : ||x|| < \\x — v \\ for all v G L}. 
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radius, we get an upper bound for the lattice sphere covering problem, which is the 
problem of minimizing the covering radius among the lattices of fixed determinant 
(see IICS991 Chapter 2] and BSV0610 . The computation of the covering radius of the 
Leech lattice in MCS99I Chapter 23] had a major impact on the study of hyperbolic 
reflection groups (see JCS99J Chapter 27]). An upper bound for the Frobenius 
number of a set of integers can be obtained from the covering radius of a suitable 
lattice (see IIFR0 5]). A recent application comes from public key cryptography; 
Micciancio MMic041 found a new connection between the average-case complex- 
ity of finding the packing radius and the worst-case complexity of determining the 
covering radius. In information theory, the quality of a lattice as a vector quantizer 
is determined by its quantizer constant (see HGG92[ IELZ051 ISB03I and RCS991 
Chapter 2.3, Chapter 21]). 

The structure of this paper is as follows. In Section [2] we discuss the com- 
putational complexity of the covering radius problem. We prove that the related 
problem of counting vertices of the Voronoi cell is # P-hard. As a byproduct of 
our construction, we show that the lattice isomorphism problem is at least as dif- 
ficult as the graph isomorphism problem. We turn to practical algorithms for the 
covering radius problem in Section [3] There we describe an algorithm which com- 
putes the vertices of the Voronoi cell of a lattice. Based on this algorithm we give 
an algorithm for computing the quantizer constant in Section [4] In Section [5] we 
report on computations with our implementation. We determine the exact covering 
radius and quantizer constants of many prominent lattices which were not known 
before. 

2. Computational complexity 

We formulate the covering radius problem as a decision problem. 

Problem 1. Covering radius problem 
Input: m, n, B £ Q mxn , p, £ <Q>. 
Output: Yes, if ii(B7L n ) < p,, No otherwise. 

It is conjectured (see MMic041 Section 1.1]) that the covering radius problem is 
NP-hard. Haviv and Regev MHR061 showed that there is a constant c v so that the 
covering radius in the / p -norm is Ll^-hard to approximate within a constant less 
than c p for any large enough p. In [GMR05] Guruswami, Micciancio and Regev 
proved that approximating it within a factor of 0(y / m/logm) for a lattice of rank 
m cannot be NP-hard unless the polynomial hierarchy collapses. 

Currently, there is only one known general and practical method to compute 
p{L) for a lattice L: First one enumerates the vertices of V(L) and then one finds 
the vertex with largest norm. The number of vertices of V(L) can be as large as 
(m + 1)! and furthermore, as we show in Theorem [T] even computing this number 
is # P-hard. 

Problem 2. Vertices of a lattice Voronoi cell 

Input: m,n,Be Q mxn . 

Output: Number of vertices of \ {BIT 1 ). 
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Theorem 1. The problem "Vertices of a lattice Voronoi cell" is # P-hard. 

It will be obvious from the proof that we could restrict the problem to the case 
m = n. We reduce the problem "Acyclic orientations of a graph", which Linial 
HLin861 showed to be # P-complete, to Problem |2j 

Problem 3. Acyclic orientations of a graph 
Input: A graph G = (V, E). 

Output: The number of orientations of G with no directed circuit. 

The structure of the proof of Theorem Q] is as follows: In Section |2~T1 we con- 
struct a matrix B with columns indexed by E defining a lattice L{G) = BZ E 
from G in polynomial time. Then we show that the vertices of the Voronoi cell 
of V(L(G)) are in bijection with the acyclic orientations of G. To establish this 
bijection we need several intermediate steps. In Section |Z21 we associate to G a 
hyperplane arrangement H{G) whose chambers are in bijection with the acyclic 
orientations of G. In Section 12.31 we recall that the chambers of a hyperplane ar- 
rangement are in bijection with the vertices of a zonotope associated to the hyper- 
plane arrangement. These two steps are standard and we cover them rather briefly. 
In Section [2~4l we show that the Voronoi cell of L(G) is a zonotope which, up to a 
linear transformation, is the one associated to the hyperplane arrangement H{G). 
In Section 12.51 as a byproduct of this construction, we show that the lattice iso- 
morphism problem is at least as difficult as the graph isomorphism problem. Some 
related complexity results concerning vertex enumeration of polyhedra given by 
linear inequalities are in MKBBEG081 |Dy83[ . 

2.1. From graphs to lattices. Let G = (V, E) be a connected graph with vertex 
set V = {1, . . . , n} and edge set E. We consider the following orientation of the 
edges of G: The head of an edge e = {v , w} G E is e + = max{f , w} and the tail 
is e~ = mm{v, w}. 

Let T C E be the edge set of a spanning tree of G, and let e G T. Deleting 
e from T divides T into two connected components with vertex sets T e + and T~, 
where e + G T+ and e~ G T~. Define the vector 6y,e G % E by 

' 1, if /+ G r+and/" G T~, 

b T , e (f) = { "!> if r e T+ and /+ G T~, 
0, otherwise. 

V 1 



Then 



L(G,T) = |^a e 6 Tie :a e Gzl C 

I eeT J 



is a lattice of rank n — 1. 

Proposition 1. Let T and T' be spanning trees ofG. Then, L(G, T) = L(G, T'). 

Proof. Since one can connect any two spanning trees by a sequence of transfor- 
mations of the form T<-^T\{e}U{/}it suffices to prove the proposition for 
T' = T\{e}U{/}. Let g G T'. If g = f, then b T ,j = ±b T , e - HgST, then denote 
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by C the cycle containing e and f.lfg^C then br> >g = &T,g- The subgraph of G 
with edge set T \ {e, g} has three connected components, denoted by C\, C%, C3. 
Given h = {v,w} 6 E, the value of br', g (h), bT, g {h) and bT, e {h) depends only 
on which connected component v and w belong to. So, in computing br', g , we can 
reduce ourselves to the case when G is the complete graph on {1, 2, 3}, g = {1, 3}, 
e = {1, 2} and / = {2, 3}. Easy computation gives bT', g = br, g + &T,e and so we 
conclude that br\ g = bT, g + abx, e with a G {— 1, 0, +1}. □ 

In the following we omit the spanning tree T from the notation L(G, T) and just 
write L(G). Note that one can find a basis of L(G) given G in polynomial time. 

2.2. From graphs to hyperplane arrangements. A matrix 

(l) V=(v 1 ,...,v m )eR nxm 

with non-zero column vectors i>j £ ]R n gives an arrangement of hyperplanes 

H{V) = {H U ..., H m } with Hi = {c€R n :c-Vi = 0}. 

The hyperplane arrangement H(V) divides the space R n into polyhedral cones, 
called regions, of different dimensions. The regions are partially ordered by inclu- 
sion and full-dimensional regions are called chambers. 

To associate a hyperplane arrangement H(G) with G we consider the incidence 
matrix Dq G M v ' xE of G which is given by 

( 1, if« = e+ 
D G (v,e) = < -1, if v = e~, 
[ 0, otherwise. 

Then we define the hyperplane arrangement of G by H(G) = H{Dq). 

In MGZ831 Lemma 7.1] Greene and Zaslavsky show that the chambers of H(G) 
are in bijection with the acyclic orientations of G: Let E be an acyclic orientation 
of E. Then a chamber of H{G) is given by 

Reg(E) = {x G R v : x v < x w if (v, w) G E}. 

Let R be a chamber of H(G). Then an acyclic orientation of E is given by 

E(R) = {(v, w) : {v, w} G E and < x w for every x G R}. 

Obviously, Reg(E(R)) = R. 

2.3. Hyperplane arrangements and zonotopes. The matrix V in (OQ) defines a 
zonotope Z(V) by 

^(V r ) = |^a i ^:-l<o i <l|. 

The faces of Z(V) are partially ordered by inclusion. It is a well-known fact (see 
e.g. BZie951 Theorem 7.16]) that the partially ordered set of regions of the hyper- 
plane arrangement H(V) is anti-isomorphic to the partially ordered set of faces 
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of Z(V): Let R be a region of H(V). Let x G R. Then the corresponding face 
Eace(.R) of Z(V) given by 



Face(i?) = < y G Z(V) : x ■ y = max x ■ 
{ z&Z(V) 

does not depend on the choice of x. Let F be a face of Z(V). Let y be in the 
relative interior of F. Then the corresponding region Reg(F) of H(V) given by 



Reg(-F) = < x G K n : max x ■ z = x ■ y 

does not depend on the choice of y. Obviously, Face(Reg(F)) = F and F 1 Q F 
if and only if Reg(F') D Reg(F). In particular, the chambers of -H"(V) are in 
bijection with the vertices of Z(V). 

2.4. From lattices to zonotopes. Let L C IR n be a lattice. The support of a vector 
w G L is v = {2 G {1, . . . , n} : 7^ 0}. The vector i> is called elementary if 
u € {— 1, 0, +1}™ \ {0} and if v has minimal support among all vectors in L \ {0}. 
We say that two vectors v, w G L are conformal if UjU), > for alH = 1, . . . , n. 
The lattice L is called regular if for every vector u G L \ {0} there exists an 
elementary vector u G L with m C u. 

Lemma 1. ( HTut711 Chapter 1]) 

(i) For any graph G the lattice L{G) is regular. 

( ii) IfL is a regular lattice, then every v G L can be written as a sum ofpairwise 
conformal elementary vectors. 

( Hi) If L is a regular lattice, v G L is elementary, and u G L satisfies u = Vj 
then there exists a factor a G Z such that u = av. 

A vector v G L for which V(L) n {x G W 1 : x ■ v = \v ■ v} is a facet of V(L) 
is called relevant. Voronoi characterizes in HVor081 page 277] the relevant vectors 
of L: A nonzero vector v G L is relevant if and only if ±v are the only shortest 
vectors in v + 2L. 

Proposition 2. In a regular lattice, a vector is elementary if and only if it is rele- 
vant. 

Proof. Let v G L be a relevant vector. By Lemma Q] (ii), we can write v = 
X^fcLi as a sum of pairwise conformal elementary vectors Wk G L. Assume that 
m > 2. Defining u = v — 2w\ gives n 7^ ±t> and u-u = v-v—A{v—W\)-wi. Since 
the vectors u^, A; = 1, . . . , m, are pairwise conformal we have (v — w±) ■ w\ > 0, 
and ±v is not the unique shortest vector in v + 2L. In this case v cannot be a 
relevant vector. Hence, m = 1 and v is an elementary vector. 

Let v G L be an elementary vector, and let u G v + 2L be a lattice vector with 
it ^ ±v. We have u — u G 2L C 2Z n and G {—1,0, +1}, which shows v C u. 
The case v/u immediately leads to f • u < n • u. If u = n, then by Lemma [U(iii), 
there exists a factor a E Z\{-1,+1} so that u = av, hence v ■ v < u ■ u. In both 
cases ±v are the only shortest vectors in v + 2L. Hence, v is a relevant vector. □ 
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The following special case of the Farkas lemma is proved e.g. in MRoc70l Theo- 
rem 22.6]. 

Lemma 2. Let L C W 1 be a regular lattice. Let x G M n be a vector, and let 
a>i, . . . ,a n G iU {±oo}. Either there exists a vector y' G (linL)- 1 - lying in 
x + nr=it — a *' a i]> or there exists a vector y G linL such that for all z G x + 
nr=it — a *' a i\ tne inequality y ■ z > holds. If the second condition holds, then 
one can choose y to be an elementary vector of L. 

Theorem 2. Let L C W 1 be a regular lattice. Let P G ~R nxn be the matrix of the 
orthogonal projection of «" onto linL. Then, V(L) = \Z{P) = P([-l/2, l/2] n ). 

Proof. Suppose that x G [-1/2, l/2] n . For all v G IT 1 \ {0} the inequality x ■ v < 
if • v holds. Write x = y + y' with y = Px G linL and y' G (linL)- 1 . For all 
wGL \ {0} we have y ■ v = x ■ v — y' ■ v < ■ v. Thus, Px G V(L). 

Suppose now that y G V(L). If there exists x G (-y + [-l/2, l/2] n )n(linL)- L , 
then y + x G [—1/2, l/2] n and P(y + x) = y. Assume that such a vector does 
not exist. Then by Lemma |2] there is an elementary lattice vector v G L so that 
v- {~y+ [-1/2, l/2] n ) > 0. This implies v ■ (-y - \v) > 0. Hence, -y V(L). 
Since V(L) is centrally symmetric, this contradicts the assumption y G V(L). □ 

In |Big97| Proposition 8.1] Biggs shows that for the lattice L(G) the matrix P 
can be written in the form P = XDq where Dq G My xE is the incidence matrix 
of G and X G ~R. ExV is given by 

number of spanning trees T with e G T and v G T+ 

(2) X(e, v) = B : — £-. 

number of spanmng trees of G 

Furthermore, the linear map given by X restricted to the image of Dq is a bijection. 

Thus, the zonotope Z(P) which is the Voronoi cell of L{G) equals ^XZ(Dg)- 
Hence, there is a linear isomorphism between the faces of V(L(G)) and those of 
Z(Dq). This completes the proof of Theorem [TJ 

Using a straightforward computation we get the following proposition. 

Proposition 3. Using the notation in $2$, the covering radius of the lattice L{G) 
is given by 



(3) ML(G)) 2 = 



-vav^j e< - E y f€E j 



Unfortunately, we do not have a combinatorial interpretation of ©. Finding one 
could lead to a proof of the NP-hardness of the covering radius problem. 

2.5. Lattice isomorphism problem. Using the construction L(G) used in the 
proof of Theorem [T] we reduce the graph isomorphism problem to the lattice iso- 
morphism problem in polynomial time. We don't know whether one can give a 
reverse polynomial time reduction. For the graph isomorphism problem no poly- 
nomial time algorithm is known. It is generally believed to lie in NP n co-NP. 
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So it is unlikely that it is NP-hard. For more information on the computational 
complexity of this problem, see the book MKST931 of Kobler, Scheming and Toran. 

Problem 4. Lattice isomorphism problem 

Input: m, re, B, B' € Q mxn matrices of rank m. 

Output: Yes, if there is an orthogonal transformation O so that OBJ/ 1 = 
B"L n , No otherwise. 

Problem 5. Graph isomorphism problem 

Input: Graphs G = (V, E G ), H = (V, E H ). 

Output: Yes, if there is a permutation a : V — > V so that for all v,w € V 
we have {v, w} 6 Eq if and only if{a(v), cr(w)} € Ejj, No otherwise. 

Theorem 3. There is a polynomial time reduction of the graph isomorphism prob- 
lem to the lattice isomorphism problem. 

Proof. Let G = (V, E G ) and H = (V, E H ) be graphs. We modify G and H by 
adding three extra vertices to V each adjacent to all vertices in V. We call the new 
graphs G' and H' which are by construction 3-connected and they are isomorphic 
if and only if G and H are isomorphic. 

It is clear that the lattices L(G') and L(H') defined in Subsection 12. II are iso- 
morphic whenever G' and H' are. For this direction it would be enough to consider 
the original graphs G and H. 

Now suppose that the lattices L(G') and L(H') are isomorphic. We apply the 

2- Isomorphism-Theorem of Whitney (actually we only use the easy subcase of 

3- connected graphs 110x1921 Lemma 5.3.2]): Because the graphs G' and H' are 
3-connected and there is a bijection between the elementary vectors preserving 
conformality, the graphs G' and H' are isomorphic. □ 



In this section we describe an algorithm which computes all vertices of a lattice 
Voronoi cell. Our focus is on implementability and practical performance, using 
the symmetries of the lattice. In fact, the algorithm computes all full-dimensional 
Delone cells and the adjacencies between them up to equivalence. We give nec- 
essary definitions in Section [3TT1 In Section [3721 we describe the algorithm's main 
steps and in the following sections we give details about its subalgorithms. In 
Section 13.61 we explain how to use Gram matrices instead of lattice basis and in 
Section [3771 we compare our method with existing algorithms. 

3.1. Notation. From now on, we assume lattices L C W l to have full rank n. 

To encode the vertices of V(L) we use Delone cells. A point x € W 1 defines a 
Delone cell D(x) by 



Denote by S(x, r) the sphere with center x and radius r. For r = min^gx, ||x — u||, 
the sphere S(x, r) is called empty, since there is no lattice point inside. In this case 



3. Algorithms 



D(x) = conv < v G L : \\x — v\\ = min \\x — w 
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the polytope D(x) is the convex hull of S(x, r) n L. The Delone cell of a vertex 
of V(L) is characterized among all Delone cells by the following properties: The 
origin is a vertex of D(x) and D{x) is full-dimensional. 

It is well known (see e.g. IIEdeOllO that the Delone cells are the projections of 
the faces of the infinite (n + 1) -dimensional polyhedral set 

Lift(L) = conv{(x, ||a;|| 2 ) : x € L} . 

The task of finding a vertex of a Delone cell of a point x, given a lattice basis 
of L, is called the closest vector problem. Generally this is an NP-hard problem 
HDKS031 ; however, there are algorithms and implementations available which can 
solve this problem rather fast in low dimensions. 

The orthogonal group O(L) of L is the group of all orthogonal transformations 
A € 0(R n ) fixing L, i.e. A(L) = L. The isometry group Iso(L) of L is the group 
generated by O(L) and all lattice translations t v : W 1 — > M. n with t v (x) = x + v 
for v € L. 

We say that two vertices x and x' of V(L) are equivalent if there is an A £ O(L) 
so that A{x) = x' . Correspondingly, we say that two Delone cells D(x) and D(x') 
are equivalent if there is an A E Iso(L) so that A(D(x)) = D(x'). 

3.2. Main algorithm. Our algorithm finds a complete list of inequivalent full- 
dimensional Delone cells of L with respect to Iso(L). The enumeration process is 
a graph traversal algorithm of the graph of equivalence classes of full-dimensional 
Delone cells of L. Two equivalence classes are adjacent whenever there is a facet 
between two of its representatives. Note that this graph can have loops and multiple 
edges. 

For the graph traversal algorithm below one needs four subalgorithms, which 
we explain in the following sections. 

Input: n, B G Q nxn matrix of rank n. 

Output: Set M. of all inequivalent full-dimensional Delone cells of the lattice 
B7L n with respect to the group lso(BZ n ). 

x <— an initial vertex of \{B'L n ). (Section |33T ) 

T <— {D(x)}. 

while there is a D E T do 
M ^MU{D}. 
T ^T\{D}. 

T <— facets of D. (Section [ 

for F 6 T do 

D' <— full-dimensional Delone cell with F = D PI D' . (Section [ 

if D 1 is not equivalent to a Delone cell in M. U T then (Section [ 

T<-TU{D'}. 
end if 
end for 
end while 
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Two full-dimensional Delone cells D(x) and v + D(x), both containing the 
origin, are equivalent under O(L) if and only if and — v are equivalent under the 
stabilizer group of D(x) in Iso(L). As a consequence, we can compute the vertices 
of V(L) under O(L) in the following way: For every orbit of full-dimensional 
Delone cells given by a representative D(x), we compute the orbits of vertices 
of D(x) under the stabilizer group and get the corresponding orbits of vertices of 
V(L) under O(L). 

3.3. Finding an initial vertex. Now we explain a method for computing an initial 
vertex of the Voronoi cell of a lattice, i.e. a full-dimensional Delone cell containing 
the origin. The method we propose is a so-called cutting-plane algorithm, which is 
a well-known technique in combinatorial optimization. 

Let us describe the geometric idea. We start with an outer approximation of 
the Voronoi cell given by linear inequalities. The first outer approximation is the 
polytope defined by the inequalities ±6j • x < \bi ■ b{ for given lattice basis vectors 
b\,...,b n . Then we find a vertex x of the approximation by linear programming 
(see e.g. BSch8610 . Deciding whether the vertex x belongs to the Voronoi cell V(L) 
can be done as follows: Compute the vertices of the Delone cell D{x). If the origin 
is a vertex of D(x), then a; is a vertex of V(L). Otherwise x is not contained in 
V(L), and for all vertices v of D{x) we have the strict inequality \\x — v\\ < \\v\\. 
So the new linear inequalities v ■ x < \v ■ v together with the old ones provide a 
tighter outer approximation of the Voronoi cell. Since we started with a compact 
outer approximation, finitely many iterations of these steps suffice to find a vertex 
of the Voronoi cell. 

One advantage of this method is that the computation of all facets of the Voronoi 
cell is not required, i.e. we do not use Voronoi's characterization (see Section I2l4l) 
of facet defining vectors, which involves solving exponentially many closest vector 
problems. Figured] illustrates this algorithm. 

Input: n, B = (bi, . . . , b n ) £ Q nxn matrix of rank n. 
Output: vertex x of Y{BZ n ). 

c <— random vector in Q n . 




Figure 1. Finding an initial vertex of V(L) 



B^{±h,...,±b n }. 
do 

x <— a vertex of the polytope {x : b ■ x < \b ■ b for all b € B}, 
which maximizes c • x. 
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£■ <— set of closest lattice vectors in BZ n to x. 
if € £ then 

return x. 
end if 
B <- B U £. 
end do 

3.4. Computing facets of, and finding adjacent Delone cells. We want to de- 
termine the facets of a full-dimensional Delone cell, which is given by its vertex 
set. This representation conversion problem can be solved by many different meth- 
ods. For details and implementations we refer to cdd HFuk95ll . Irs HAvi93ll . pd 
HMar97H and porta ICT97l . 

In order to exploit the symmetries we use the adjacency decomposition method 
(see HCR96llBDgQ7llDSV07l ). It allows to compute a complete list of inequivalent 
facet representatives: We compute an initial facet by linear programming and insert 
it into the list of orbit representatives of facets. From any such orbit, we compute 
the list of facets adjacent to a representative and insert it, if necessary, into the list 
of representatives until all orbits have been treated. Computing adjacent facets is 
itself a representation conversion problem in one dimension lower. So this method 
can be applied recursively (see MBDS071 IDSV07ID . Note that our main algorithm 
is itself an adjacency decomposition method. 

After the computation of facets, we can compute adjacent full-dimensional De- 
lone cells: We take an initial vertex v and so get a tentative empty sphere. If the 
sphere is not empty, then we find another vertex v and iterate until the sphere is 
indeed empty. Figure |2]illustrates this algorithm. 

Input: n, B € Q nxn matrix of rank n, a full-dimensional Delone cell D and a 
facet F of D. 

Output: vertex set V' of a full-dimensional Delone cell D' with D n D' = F. 

(f> <— affine function on R" with F = {x G D : <f>(x) = 0} and 4>{x) > 

on D - F. 
Vf <— vertices of D belonging to F. 
v *— a point of BZ n with <f>(v) < 0. 
do 
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S(c, r) <— sphere around Vf U {v}. 
V 1 <— closest vectors in BZ n to c. 
if — c\\ = r for au'eV' then 

return V'. 
end if 

v <— one element of V'. 
end do 

One way to speed up the convergence of this algorithm in practice is to heuristi- 
cally choose an initial vector v with a sphere S(c, r) of small radius. 

3.5. Checking equivalence. We have to test equivalence and compute stabilizers 
under the group Iso(L) of Delone cells of different dimensions. Below we propose 
three different methods for this. 

We can encode a Delone cell D by the center c(D) of the empty sphere around 
it or by the vertex barycenter g(D) = p^-pi YlvevertD v °f i ts vef tex set vert D. 
Both c(D) and g(D) are invariant under the stabilizer of D. Any two full-dimensional 
Delone cells D, D' are equal if and only if c(D) = c(D'). However, it is possible 
if n > 3 that c(D) lies outside or on the boundary of D. If c(D) lies on the bound- 
ary of D then a facet containing c(D), which is itself a Delone cell, has the same 
center as D. Hence, the sphere centers can be used to distinguish full-dimensional 
Delone cells but they do not distinguish Delone cells. Therefore, we use the vertex 
barycenter. 

In the first method we consider the classes of the vertex barycenters g(D) and 
g(D') in the quotient K n /L and check their equivalence under the finite group 
Iso(L)/£, ~ O(L). The generic methods underlying isomorphism and stabilizer 
computations generate the full orbit of g(D) under Iso(L)/^. This is typically 
memory intensive. In some cases we can use a method from computational group 
theory, which we now explain in an example. Suppose g(D) is expressed in a 
basis (pi, ... , b n ) of L as . . . , ^) with < aj < 5 and we want to compute 
its stabilizer under the group Iso(L)/£. The vector 2g(D) reduced modulo L is 
expressed as (^-, . . . , ^) with < dj < 2. We first compute the stabilizer H 
of the vector 2g(D) under the action of Iso(L)/^. The stabilizer of g(D) under 
Iso(L)//, is equal to the stabilizer of g(D) under H. This method generalizes to 
more than two prime factors and it is more memory efficient because the generated 
orbits are smaller. 

The second method uses finite metric spaces of the vertex set of full-dimensional 
Delone cells obtained from the metric \\v— v' || 2 HDL97I Chapter 14]. A finite metric 
space defines an edge-weighted graph. Testing if two edge weighted graphs are 
isomorphic can be reduced to testing if two vertex-weighted graphs are isomorphic 
(see MMcK061 Page 25]). In practice, the program nauty [McK06] can solve the 
isomorphism problem if the number of vertices of D and D 1 is not too large. If the 
metric spaces are not isomorphic, then D and D' are not equivalent under Iso(L). If 
they are isomorphic then every graph isomorphism corresponds to a linear isometry 
between D and D' MBDS07llDSV07i For each of those isomorphisms, we check if 
it belongs to Iso(L). This method is useful when the isometry group of D is small. 
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For the third method, we use laminations over the ?7,-dimensional lattice L. 
Let D be a Delone cell of L with vertex barycenter g. One defines an (n + 1)- 
dimensional lattice L(g) by embedding L C W 1 into ]R™ +1 and adding layers to 
it 



where v G M n+ is chosen so that v + g is orthogonal to the space spanned by 
L and normalized so that ||u + g\\ = 1. A variant of this construction is used for 
example to build the laminated lattices, see HCS991 Chapter 6]. If <fi is an element 
of 0(L(g)) preserving every layer of the lamination, then it maps the vector v to 
some vector w = v + h with h G L. The function x h- > (/>(x) + /i preserves the 
Delone cell and every element preserving the layers is obtained in this way. In 
practice, we can use the program AUTO (see HPP97II ) of the package CARAT (see 
HOPS981D for computing this automorphism group. The isomorphism problem is 
treated similarly using the program ISOM. 

3.6. Working with Gram matrices and periodic structures. In many cases, it is 
more convenient to work with the Gram matrix B B instead of the lattice basis B 
(see HCS99I Chapter 2.2]). For instance, when B is irrational but B B is rational. 
Note that our algorithms can be reformulated in terms of Gram matrices. Note also 
that all our algorithms can be modified to deal with periodic point sets, that is for 
finite unions of lattice translates. Our implementation is available from |Dut08l . 

3.7. Comparison. In HVB961 Viterbo and Biglier describe another algorithm for 
computing the Voronoi cell of a lattice, called the diamond cutting algorithm. As in 
our approach they start with a parallelepiped P defined by the basis vectors. Then 
they determine all lattice vectors which lie in a sphere containing 2P. This set 
contains all facet defining lattice vectors of V(L). Successively they add cutting 
planes obtained from these vectors and update the complete face lattice of the ten- 
tative Voronoi cell. They terminate when its volume coincides with det L. Their 
implementation uses floating point arithmetic. 

In comparison, our approach has the following advantages: We use the presence 
of symmetry in an efficient way. We do not need to compute a huge initial list of 
potential facet defining lattice vectors. Our algorithm does not need to compute the 
face lattice, not even for computing the quantizer constant as explained below. Our 
implementation uses rational arithmetic only. 



Recall from the introduction that the quantizer constant of a lattice L is the 
integral 



JV(L) 

A standard method for computing the integral G{L) is to decompose V(L) into 
simplices. Suppose that .Sis a simplex with vertices v\, . . . , v n+ \ in W 1 . Then (see 



L(g) = {av + h : h G L,a G Z} CM 



n+l 



4. Computing quantizer constants 
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Thus G(L) can be obtained by summing the integrals of all simplices in a decom- 
position of V(L). Several practical methods for decomposing a polytope into sim- 
plices are discussed in MBEF 981. In our implementation, we use the triangulation 
obtained by the program Irs. However, this method as well as the other methods 
explained in HBEF981 are sometimes impractical and they do not use symmetries. 

In order to get a group invariant decomposition, we can use the barycentric sub- 
division of P. That is, given any flag Fo C F\ C . . . C F n = P of faces of 
P, we associate the simplex with vertex set go,9i, ■ ■ ■ ,9n where g\ is the vertex 
barycenter of F, L . Note that in general there is a difference between the bary center 
^-p J p xdx of a polytope P and its vertex barycenter |vcr \ p| Revert p v - The 
group acts on the barycentric subdivision and the stabilizer of each simplex is triv- 
ial. In practice, the number of orbits of flags can be too large. 

We propose a hybrid approach, which combines the benefits of both methods. 
Let T be the facet set of an n-dimensional polytope P. We can assume without 
loss of generality that P has the origin as its vertex barycenter. We then have 



To compute this sum, it is sufficient to compute the integrals only for orbit repre- 
sentatives of facets. Let F be a facet of P and pp be a point in the affine space 
spanned by F. Then we can transform the integral over the cone conv(F, 0) in the 
following way: 



If Pi? is the orthogonal projection of the origin onto F then the second summand 
vanishes. This point may not be invariant under the automoiphism group of the 
facet F, but the vertex barycenter is. If we use the vertex barycenter, we also 
have to compute the barycenter of the polytope F as well as the volume and the 
square integral. In order to use symmetries coming from non-orthogonal linear 
transformations of P, we use the matrix valued integral 



0991 Chapter 21, Theorem 2]) the following holds: 




(4) 






This integral splits according to 
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where 

J (P)= / dx = volP, h{P)= / xdx, h(P)= / xx l dx. 
Jp Jp Jp 

Let G be a group of automorphisms of P. If g € G acts on R n as x i— > Ax + v then 

we define H (g) = ( . ) the corresponding (n + 1) x (n + 1) matrix acting 



^ v A f 

on homogeneous coordinates. Let 0±, . . . ,O r be the G-orbits of facets of P, with 
representatives Fx, ... , F r . Then the integral io,i,2(-P) simplifies to 

Io,i, 2 (P) = E ( T^T E^Ko,i,2(conv(F,0) 
i=i VJ 1 9 eG 

Assume that Io,i,2(conv(Fj, 0)) is already computed. To compute the sum in the 
parenthesis, we first incrementally compute a basis of the affine hull of the orbit 
{H(g)Io t i t 2(conv(Fi, 0))H(g) t : g G G}. The only G-invariant element of the 
affme hull is the sum we want to compute. 

We now want to compute Jo,i,2(conv(F, 0)) in terms of lower dimensional inte- 
grals. The integral depends on the chosen basis. If / is an affine transformation of 
W 1 , then the change of variables formula for integrals gives 

H(f)I , 1>2 (fP)H(fy\detH(f)\=I 0A>2 (P), 

for any an n-dimensional polytope P in W 1 . This allows to compute Io,i,2(-P) for 
another basis. So, we can choose a coordinate system such that 

1 : x G F' !> C 

x 

where F' C W 1 " 1 is an (n — 1) -dimensional polytope. We then have the following 
formulas: 

/ (conv(F, 0)) = i/o(F'), ^(conv^, 0)) = ^ ( 

I 2 (conv(F,0)) = ^2/o,i,2(F')- 

For computing Io,i,2(F')' we have two options: Either we use the first method 
of this section, which involves computing a triangulation or we apply the above 
method recursively. The decision is made heuristically, depending on the size of 
the automorphism group of F and its number of vertices. In order to reduce the 
size of the computation, one can store intermediate results. 

Those methods are general and apply to any polytope and any polynomial func- 
tion, which we want to integrate over P. Note that a similar method of using the 
standard formula (Tjfl) has been used for computing the volume in [B EF981 under 
the name of Lasserre's method ( BLas98fl ). albeit in a non-group setting. 
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5. Results 

In this section, we collect results from our implementation of the algorithms 
explained in Sections [3] and 01 We obtain previously unknown exact covering den- 
sities and quantizing constants of several prominent lattices and their duals. Recall 
that the dual L* of a lattice L C W 1 is defined by 

L* = {x G R n : y ■ x G Z for all y G L} . 

The covering density of an n-dimensional lattice L is 

where B n is the unit ball in R n . Other computations of Voronoi cells of lattices 
can be found in 09TI HEMS03I Chapter 5] and HMP95H . All computations are 
done in exact rational arithmetic. In the tables the covering densities are given in 
floating point; the exact expressions would be too large. 



5.1. Coxeter lattices. The root lattice A n is defined by 

A n = |xGZ n+1 : |>i = oj. 
If r divides n + the Coxeter lattice A n (see BCox51IO is defined by translates of 

An- 

A n = A n U « + A n ) U . . . U ((r - l)v r n + A re ), 

where v r n = Y^i=2( e i ~ e i)- The dual lattice of A n is J\n +1 ^ r . 

The Delone decomposition of the lattice A n has been studied in MAnz021lAnz061 
|Bar941 up to dimension n = 15, hoping to obtain lattices with low covering 
density. One pleasant fact is that the symmetry group of A n contains the group 
Sym(n + 1) x TL^. Latter can be represented as a permutation group acting on 
n + 3 points, which drastically simplifies isomorphism computations. 

In Table Q] we list the obtained results. Note that the lattices Af 7 , A}g, A| and 
A3,} turn out to give new record sphere coverings. Up to dimension 8 all those lat- 
tices are well known and their Voronoi cells can be obtained by standard computer 
algebra software. Our list is complete up to dimension 21. For the missing cases 
we could not finish the computation. 



5.2. Laminated lattices. The laminated lattices, which are defined in IICS991 Chap- 
ter 6], give the best known lattice sphere packings in many dimensions. The Delone 
subdivision is known up to dimension 8 and in dimension 24 for all laminated lat- 
tices and their duals ICS991 Chapters 21, 23, 25]. In dimension 16, the covering 
density of A 16 is known IICS991 Chapter 6]. 

In Table|2]we list the obtained results, which are complete up to dimension 17. 
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lattice 


# orbits 


covering density 


lattice 


# orbits 


covering density 




6 


18.543333 


A 9 


5 


4.340184 


A 2 

A 4 


6 
16 


94.090996 
5.598337 


A 3 

A 

A ii 


11 
4 


27.089662 
7.618558 


A 2 

M 13 


10 


134.623484 


A 7 

M 13 


10 


7.864060 


A 3 


17 


32.313517 


A 5 


31 


9.006610 


A 2 

A 8 
A 15 


10 
10 


722.452642 
11.601626 


A 4 


19 


25.363859 


A 2 
A 17 

fk\ 7 


15 

73 


1068.513081 
12.357468 


A 3 
M 17 

A 9 
M 17 


26 
24 


240.511580 
17.231927 


A 2 

A 5 
A 19 


15 

80 


5921.056764 
25.609662 


A 4 

A 10 

M 19 


58 
80 


40.445924 
21.229200 


A 3 
A 20 


40 


307.209487 


A 7 
A 20 


187 


20.366828 


A 2 
A 21 


21 


8937.567486 


A 11 
A 21 


64 


27.773140 


A 3 
A 23 

A 23 

A 12 
rt 23 


55 
187 
100 


2405.032746 
79 575330 
43.231587 


A 4 
A 23 

A23 


85 
495 


205.561225 
31 858162 


A 5 


144 


115.011591 








A 1 Q 

A J^-J 

rt 25 


210 


54.472182 








A 3 

A 26 


75 


3184.1387034 


A 9 

A 26 


1231 


50.937168 


A 4 
M 27 


156 


350.137031 


A 7 
rt 27 


650 


81.869181 


A 14 
rt 27 


338 


76.909712 








A 3 
rt 29 

A 6 
rt 29 

Al5 
M 29 


102 
711 

678 


25664.644103 
154.329831 
114.084219 


A 5 
rt 29 

A 10 
rt 29 


347 
3581 


202.040331 
84.324725 


A 16 
rt 31 


1225 


33.934941 









Table 1 . Number of orbits of full-dimensional Delone cells and 
covering density for some Coxeter lattices 



5.3. Shorter Leech lattice. The 4600 shortest vectors of A23 of define a sublattice 
of index 2, called the shorter Leech lattice 2 3 ( 1CS99I Page 179, 420, 441]). The 
Delone decomposition (see Table [3]) is remarkable in many respects: There are 
only 5 orbits and the first one has the full symmetry group of the lattice. It turns 
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lattice 


# orbits 


covering density 


lattice 


# orbits 


covering density 


A 9 


5 


9.003527 


A* 


9 


9.003527 


Aio 


7 


12.408839 


A* 


21 


9.306629 


A max 
yv ll 
A min 


11 
18 


24.781167 
24.781167 


A max* 
A ll 

A min* 
A ll 


18 
153 


19.243468 
8.170432 


A max 
A 12 
A mid 

A min 
1V 12 


5 

23 
13 


30.418954 
30.418954 
30.418954 


A max* 
A 12 

A mid* 
A 12 

a min* 
A 12 


8 

52 
78 


42.728408 
19.176309 
12.292973 


A max 
iv 13 

A mid 
A 13 

A TThtTL 

Il 13 


18 
46 
129 


60.455139 
35.931846 
60.455139 


A max* 
A 13 

A mid* 
A 13 

A 13 


57 
125 
5683 


43.214494 
19.155991 
13.724864 


Al4 


65 


98.875610 


A* 
A 14 


1392 


34.721750 


A15 


27 


202.910873 


A* 
A 15 


108 


25.642067 


Al6 


4 


96.500266 


A* 
A 16 


4 


96.500266 


An 


28 


197.719499 


A* 


720 


100.173101 


Al8 


239 


301.192334 








A 23 


709 


7609.03133 









Table 2. Number of orbits of full-dimensional Delone cells and 
covering density for some laminated lattices and their duals 



number of vertices 


size of stabilizer group 


94208 


84610842624000 


32 


1344 


24 


10200960 


24 


1320 


24 


1320 



Table 3. Orbits of full-dimensional Delone cells of O23 



out that A23 = O23 U (v + O23) where v is the center of a centrally symmetric 
Delone cell lying in the first orbit. The covering density of O23 is 15218.062669. 

5.4. Cut lattices. The cut polytope CUT n is a famous poly tope appearing in com- 
binatorial optimization (see MDL97II ). It has 2 n_1 vertices and is of dimension 
n ( n ~ 1 > _ xhe lattice generated by its vertices is called cut lattice and is denoted 
by L(CUT n ) (see HDG95H ). The polytope CUT„ is one of its full-dimensional 
Delone cells. We list our results in Table 0] 
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lattice 


dimension 


# orbits 


covering density 


L(CUT 3 ) 


3 


2 


2.09439 


L(CUT 4 ) 


6 


4 


5.16771 


L(CUT 5 ) 


10 


12 


40.80262 


L(CUT 6 ) 


15 


112 


255.4255 



Table 4. Dimensions, number of orbits of full-dimensional De- 
lone cells and covering density of some cut lattices 



lattice 



quantizer constant 



Ao 

A 9 
A 2 

A3 

A 2 
A 3 

Afa 

Kl2 

Table 5. 



151301 
2099520 

1371514291 „ 
19110297600 ~ 

2120743 , 
^5.28 13271040 

8651427563 
1 / 2T5 S "26578125000 
4568341 ^ 
64512000 " 

452059 
^5702400 
287544281699 



V4.3 10 1325839006800 
6387657954959 



/3^9 46506442752000 



0.07206 
a 0.07176 
w 0.072166 

0.072079 
0.07081 
4 0.07174 

0.070426 
0.070494 



29183629 
412776000 

797361941 



0.070700 
ps 0.070095 



\/36567561000 

Quantizer constants of some lattices 



5.5. Quantizer constants. In Table [5] we collect some new exact quantizer con- 
stants. 

According to MAE98I . the lattice D^ is conjectured to be the optimal lattice 
quantizer. Conway and Sloane approximated G(K 12 ) ( BCS99I Table 2.3]) using 
Monte-Carlo integration; our exact computation fits into their bounds. 
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